Jamming in Granular Polymers 
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We examine the jamming transition in a two-dimensional granular polymer system using com- 
pressional simulations. The jamming density decreases with increasing length of the granular 
chain due to the formation of loop structures, in excellent agreement with recent experiments. The 
jamming density can be further reduced in mixtures of granular chains and granular rings, also as 
observed in experiment. We show that the nature of the jamming in granular polymer systems 
has pronounced differences from the jamming behavior observed for polydisperse two-dimensional 
disk systems at Point J. This result provides further evidence that there is more than one type of 
jamming transition. 

PACS numbers: 45.70.-n,81.05.Rm,83.80.Fg 



Jamming, or the development of a resistance to shear, 
is a phenomenon that occurs when a disordered assem- 
bly of particles subjected to increasing density, load or 
other perturbations exhibits a transition from a liquid- 
like state that can flow to a rigid state that acts like a 
solid under compression. Tremendous recent growth in 
this field has been driven by the prospect that jamming 
may be associated with universal properties across a wide 
class of systems including granular media, foams, emul- 
sions, colloids, and glass forming materials Q. One of the 
most accessible routes for exploring the jamming transi- 
tion is gradually increasing the density of a sample in the 
absence of shear or temperature. Here, jamming occurs 
at a density termed 'Point J.' Jamming transitions 
have been studied both for noncohesive granular media 
and for cohesive and/or nonspherical granular materials 
There is considerable evidence that for frictionless 
disordered disk assemblies, critical behavior occurs near 
Point J, with both the pressure and the particle coor- 
dination number Z exhibiting power law behavior as a 
function of packing density </) 0, 043 ■ Similar behavior 
appears when the shear, external forcing, or temperature 
are finite, providing further evidence that the jam ming 
transition may indeed exhibit universal properties |10l - 
16[. If such universal behavior holds for other systems 
that undergo jamming, it would have profound implica- 
tions for the understanding and control of disordered and 
glassy systems. 

The most widely studied two-dimensional (2D) jam- 
ming system contains bidisperse frictionless disks. When 
two sizes of disks with a radius ratio of 1 : 1.4 are mixed 
in a 50:50 ratio, a jammin g tra nsition occurs at a density 
of </> = 0.84 0, MJM- To explore whether the 



jamming transition is universal in nature, it would be 
ideal to have a system in which the jamming density 0c 
could be tuned easily. Here we propose that one model 
system which meets this criterion is an assembly of 2D 
granular polymers. This model is motivated by experi- 
ments on granular polymers or chains of the type used 
for lamp pulls, where various aspects of knot formation. 



diffusion processes, and pattern formation have been ex- 



plored |17Hl9j . We model the chains as coupled har- 



monically repulsive disks similar to those studied in the 
polydisperse disk system, with a constraint on the mini- 
mum angle that can be spanned by a string of three disks. 
Other workers have considered freely-jointed chains f^Sl 
or chains of sticky spheres 2l| . Although our model is 2D 
and neglects friction, we show that it captures the same 
features found in recent three-dimension al ( 3D) granu- 
lar polymer compaction experiments 22|, [23| . To study 
the jamming transition, we construct pressure versus (f> 
curves by compressing the chains between two walls, and 
compare our results to compression experiments on 2D 
polydisperse disks 0, [l^l- We show that the jamming 
density (j)c decreases with increasing chain length and 
saturates at long chain lengths, in agreement with the 



experiments of Ref. [22| . The decrease of (/>c occurs due 



to the formation of rigid loops along the chains which sta- 
bilize voids inside the packing. Unlike the bidisperse disk 
system where the pressure scales nearly linearly with (/>, 
in the chain system the pressure increases with a power 
law form with an exponent significantly larger than 1 as 
the jamming transition is approached. Jamming of our 
frictionless granular chains shares several features with 
jamming of frictional disks and could be distinct from 
the jamming transition for frictionless disks. 

We simulate a 2D system confined by two walls at 
X — and x — L and with periodic boundaries in the 
y direction. The wall at cc = L is held fixed while the 
position of the other wall is allowed to vary in order to 
change the density. The system contains N chains or 
loops, each of which is composed of Afb individual disks 
that are strung together by harmonic springs and that 
experience a constraining force which limits the bending 
radius of the chain. In loops, the two ends of a chain are 
connected together. The disk-disk interaction is modeled 
as a stiff harmonic repulsion, and the motion of all disks 
is taken to be overdamped in order to represent the fric- 
tional force between the disks and the underlying floor. 
A given disk i moves according to the following equation 



Here we take the damping constant rj — 1. The disk- 
disk interaction potential is F^^ — J^f^^i'' ^^(^off 

^ij I 
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)fy where the spring constant kg — 300, 
(Ri — Rj )/i?,y , Q is the Heaviside 



step function, and r^s — 



where 



is the radius 



of disk For the chains and loops we set r,; = 1; for 
a bidisperse disk system we set = 1 for half of the 
disks and n — 1.4 for the other half of the disks. The 
chain interaction potential is F* — kg{rcs — rik)rik, 
and it acts only between a disk and its immediate neigh- 
bors along the loop or chain. The bending constraint 
potential F^^ = Y.i (''stiff - ra)Q{rst\s - rii)va acts be- 
tween disks separated by one chain element, with Tstiff = 
2rcff sin(6's/2) and Og — 0.827r unless otherwise noted. 
Smaller values of Og produce more bendable chains. The 
disk-wall interaction force F^ is computed by placing a 
virtual disk at a position reflected across the wall from 
the actual disk, and finding the resulting disk-disk force. 
To initialize the system, we place the chains, rings, or 
individual disks in random non-overlapping positions to 
form a low density unjammed phase, such as in Fig.[IJa). 
The X = L wall is held fixed while the other wall is grad- 
ually moved from a; = to larger x in small increments of 
5x. The waiting time between increments is taken long 
enough so that the system has sufficient time to relax to 
a state where the velocities of all particles are indistin- 
guishable from zero. 

We identify the jamming transition by measuring the 
total force exerted on the fixed wall by the packing, 
P = X^r*^*" ■ ^"^^ the average contact number 
Z — (NMh)^^ ^ ^ function of the total den- 

sity of the system defined by the spacing between the 
two walls. To determine the contact number Zi of an in- 
dividual grain in a chain, we first count the immediate 
neighbors of the grain along the chain, and then add any 
other grains that are in contact with the individual grain. 
The force P is proportional to the Pxx component of the 
pressure tensor. At the jamming transition, the pressure 
in the system becomes finite 0, while below jam- 

ming P = 0. Previous simulations on 2D disordered disk 
packings have revealed the onset of a finite pressure near 
(j)c = 0.84 which grows as P cx {(j) — 4>c)^ with ^ = a/ — 1, 
where a/ is the exponent of the interparticle interaction 
potential 0, 0] • Recent work which includes careful cor- 
rections to scaling indicates that the exponent ip — l.l 
2J|. Theoretical work on the jamming of 2D disks also 
predicts a power law scaling of the pressure versus den- 
sity Q, and indicates that the contact number Z should 
scale as Z (X {(t) — 0c)'^, with P = 0.5. Experiments using 
a combination of shear and compression on the same disk 
system found that after performing cycling to reduce the 
effect of friction, the pressure and Z both scale with the 




FIG. 1: (Color online) Granular configurations in a portion 
of the sample. The x direction runs vertically and the fixed 
wall is at the top of each panel, (a) An unjammed system 
of = 67 granular polymers with length Mi, = 16. (b) 
The jammed configuration for the same system contains voids 
which appear when the chains form ring structures, (c) The 
jammed configuration for N = 67 loops of length Mb — 16. 
The jamming density is lower than for systems of chains or 
individual disks, (d) The jammed configuration a.t (p — 0.84 
for a sample of A'^ = 1500 bidisperse frictionless disks contains 
no significant voids. 




FIG. 2: (a) The pressure P vs <j!> for a bidisperse disk system. 
Above the jamming transition at ~ 0.84, P increases nearly 
linearly with (j). (b) P vs for the granular polymer system 
with chains of length Mj, = 6 (T), 8 (A), 10 (♦), 16 (□), 
and 24 (•). As Mb increases, the onset of a finite value of P 
indicating jamming drops to lower values of (ji and the near 
linear scaling of P with is lost. 



density as power laws with ip = 1.1 and [3 = 0.495 [9|. 

We first test our simulation geometry using the bidis- 
perse individual disk system. A configuration of = 
1500 disks appears in Fig. [TJd) just above the onset of 
a finite pressure P at w 0.84. In Fig. [Ifa) we show 
that for the disk system above jamming, P increases 
nearly linearly with cj), consistent with a scaling expo- 
nent ip = 1.1 ±0.1. This indicates that our compressional 
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FIG. 3: (a) Scaling of pressure vs density, PM^'^ vs (j) — 4>c 
close to the jamming transition for chains of length Alt = 
6 (▼), 8 (A), 10 (♦), 16 (□), and 24 (•). (b) Scaling of 
{Z - Zc)Mb YS for chains of length Mb = 6, 8, 10, 16, 

and 24, with the same symbols as in panel (a). 





FIG. 4: (a) Contact number Z vs cj> for chains of length 
Mb ^ 6 (▼), 8 (A), 10 (♦), 16 (□), and 24 (•). (b) Bulk 
pressure tensors p^x (O) and Pyy (♦) vs for Mb = 16. 



geometry captures the jamming behavior found in other 
studies of bidisperse disks. 

We use the same compression protocol to study the 
jamming behavior of granular polymers, as illustrated in 
Fig. [Ua,b) for a system with TV = 67 chains that are 
each of length Mb = 16. In Fig. [UJb) we plot P versus 
(j) for granular polymer chains of lengths Ah = 6, 8, 10, 
16, and 24. For the chains, the onset of finite P indicat- 
ing the beginning of jammed behavior occurs at a much 
lower density than for the bidisperse disk system shown 
in Fig.[2Ia), and as Mf, increases, the jamming transition 
shifts to even lower (j) and the nearly linear dependence 
of P on is lost. We illustrate scaling of P near the 
jamming transition in Fig. [SJa) where we plot PM^-^ vs 
<j) — (pc- Here we find 3. The jammed state devel- 
ops isotropic rigidity, as indicated by the plot of the bulk 
pressure tensor components Pxx and Pyy in Fig. |4ljb). To 
test whether the packing also develops a finite response to 
shear at the jamming transition, we fix the packing den- 
sity and apply a shear to the system by applying a force 
Fs/ieor = 0.04y to all particles that are in contact with 
the mobile wall. We measure the resulting steady state 
velocity Vshear — dRi/dt ■ y of all particles that are in 
contact with the stationary wall on the other side of the 
sample, discarding any brief initial transient responses. 
In Fig. [S] we plot Vshear and P versus (j) for samples with 
chains of length Adf, = 6, 8, 10, 16, and 24. In each case. 



FIG. 5: (Color online) The shear velocity Vshear of grains 
adjacent to the stationary wall vs </!> (dark black symbols) and 
the corresponding pressure P in the packing vs </> (light red 
symbols) for a sample in which a shear force is applied using 
the mobile wall with chain lengths Mb = 6 (T), 8 (A), 10 (♦), 
16 (□), and 24 (•). The onset of a finite pressure and a finite 
shear response occur at the same value of for each sample. 
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FIG. 6: (Color online) (a) The jamming threshold (pc versus 
chain length Mb for chains with a bending angle of 6s = 0.827r 
(circles) and 6s = 0.7567r (squares). (/!>c decreases with increas- 
ing Alb and saturates at large Mb- Inset: (j>c vs 6s for a system 
with Alb = 16. (b) (j>c vs chain fraction Nc/N for a system 
with Alb ~ 16 and a mixture of loops and chains. As the 
fraction of chains decreases, (he decreases. 



a finite shear response and a finite pressure P appear si- 
multaneously at the jamming density. Our use of a fixed 



P. 



at all densities rather than a variable shear rate 



affects the shape of the Vshear versus curves; variations 
in the shear rate as well as related hysteretic effects under 
shear will be the subject of a future publication. 

We define the jamming threshold 4>c as the density at 
which P rises to a finite level. The same threshold also 
appears as the sudden onset of an increase in Z, as seen 
in Fig. Ufa). The scaling of {Z — Zc)Mb vs <j> ~ (f>c ap- 
pears in Fig. IHb) , where the exponent f3 falls in the range 
(3 — 0.6 to 0.8. In Fig. [6] we plot (j)c versus Mb for chains 
of two different stiffnesses: 6s = 0.827r and 9s — 0.7567r. 
In both cases, (jjc decreases monotonically with increas- 
ing Mb and saturates for large Mb- Recent experiments 
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on the packing of granular polymers showed the same 
behavior: the final packing density decreased for increas- 
ing chain length and saturated for very long chains (2^ . 
This was attributed to the formation of rigid semiloops 
which stabilized voids in the packing and decreased the 
jamming density. Loops have also been observed in dense 
3D packings of freely-jointed chains [lH . Since the min- 
imum area spanned by a semiloop increases with 9s, the 
jamming density should be lower for larger 6^ when larger 
voids are stabilized. Figure[Hb) illustrates the voids that 
appear in our chain packings due to the formation of rigid 
semiloops. For comparison, the bidisperse disk system 
shown in Fig. [ijd) contains no large voids. If the rattler 
disks in Fig.[ljd) were removed, the amount of void space 
present would increase, but the chain system would still 
be able to stabilize a larger amount of void space since 
the constraint of the chain backbone permits the forma- 
tion of larger arches around the voids than the arches 
that can be stabilized in the bidisperse disk system. The 
semiloops increase in size for increasing Og and (pc is re- 
duced at higher 6's, as shown in Fig.jBJa). The plot of (f)c 
versus 9s in the inset of Fig.lHJa) for fixed Mi, = 16 shows 
that 0c monotonically decreases with increasing 9s. For 
perfectly flexible chains with 9s — 0, we find (t>c ~ 0.8. 
This is lower than the density of a triangular lattice due 
to the trapping of voids within the packing by the phys- 
ical constraints of the bonding between chain elements. 
If the system were annealed or shaken for a sufficiently 
long time, these voids could eventually be freed and the 
perfectly flexible chains would form a perfect triangular 
lattice. 

As a confirmation of the idea that the formation of 
rigid semiloops is the mechanism by which the jamming 
density is depressed, Ref. [22] includes experiments per- 
formed on mixtures of granular polymers and granular 
loops with equal length M5. In this case, (/)c decreased lin- 
early as the fraction of loops increased. We find the same 
effect in our 2D system by varying the number of loops 
Ni and chains Nc in a sample with fixed N = Ni + Nc and 
fixed A/fc. In Fig. [6t^b) we plot 4>c versus Nc/N, where 
Nc/N — 1 indicates a sample containing only chains and 
Nc/N = is a sample containing only loops. As Nc/N 
decreases, (f>c decreases. The jammed configuration for a 
sample with Mf, = 16 containing only loops, Nc/N = 0, 
appears in Fig.[TJc). The number of voids present is much 
larger than in the Nc/N = 1 sample shown in Fig.[IJb). 
Interestingly, the voids began to form a disordered trian- 
gular packing. 

Our results indicate that the experimentally observed 
dependence of the jam ming density on chain length or 
fraction of loops in Ref. [22| is not caused by friction or 
other possible spurious effects, but is instead a product 
of the geometrical configuration of the chains and loops. 
The surprisingly good agreement between our 2D simula- 
tions and the 3D experiments may be due to the fact that 
in each case, the semiloops formed by the chains are 2D 




FIG. 7: (Color online) (a) P vs for a chain sample with 
Mb = 20 on a log-linear scale. Light red line: initial com- 
pression; heavy black line: steady state reached after four 
compressions, (b) P vs for a system with Mi, = 38 for the 
initial compression (light red line) and for the fourth com- 
pression (heavy black line). The response is hysteretic but 
the curves remain nonlinear on all compressions. Inset: A 
portion of the P vs (j) curve during the second compression 
cycle in the same sample showing a sudden pressure drop as- 
sociated with the collapse of an unstable void. 



in nature. Additionally, in the experiment the container 
used to hold the sample induced ordering of the chains 
and loops near the walls and may have caused the system 
to act more two-dimensional. The fact that much of the 
physics observed for the 3D system can be captured in 2D 
models means that 2D experiments, which are much eas- 
ier to image than 3D experiments, could provide many of 
the same insights for understanding jamming in a system 
where the jamming density can be tuned easily. 

The chain system exhibits a pronounced hysteresis ef- 
fect that can be seen by cycling the mobile wall in and 
out to repeatedly compress and uncompress the packing. 
Figs. EJa) and (b) show a comparison of the responses 
during the first compression and during the fourth com- 
pression. After four compressions the system does not 
exhibit further hysteresis. In contrast, we find little or no 
hysteresis for the bidisperse disk system. During the ini- 
tial compression, the chain systems often exhibit sizable 
fiuctuations in P above the onset of jamming. Sudden 
drops in the pressure, such as that shown in the inset of 
Fig.EJb), occur due to the collapse of semiloops that are 
larger than the minimum stable size. After all semiloops 
have reached a stable size, we find no further hysteresis. 
Even after cycling to a steady state, the P vs curves re- 
main power law in nature with an exponent significantly 
larger than 1 and do not become linear or nearly linear. 
Simulations of compressed 2D frictional bidisperse disk 
systems show that void structures can form during the 
initial compression but collapse during subsequent cy- 
cles, allowing the sample to reach the same density as 
a frictionless disk sample In the chain system, the 
void structures are associated with semiloops that have 
formed in the chains and, unlike in the disk system, the 
voids can never be fully collapsed by repeated cycling. 

The fact that the granular polymers do not exhibit the 
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same behavior at the jamming transition as the bidis- 
perse disk systems do at Point J provides additional ev- 
idence that jamming does not occur with universal fea- 
tures in all systems, and that the criticality found in the 
bidisperse disk systems may be associated with a special 
type of jamming. Additional studies on a variety of dif- 
ferent types of systems would need to be performed to 
confirm whether the jamming behavior is indeed differ- 
ent for each system or whether there is a small number 
of different classes of jamming behaviors, with the granu- 
lar polymer system and the bidisperse disk system falling 
into separate classes. 

In summary, we have introduced a numerical model 
of 2D granular polymers that can be used to study the 
jamming transition. The onset of jamming occurs at a 
density that decreases with increasing chain length and 
saturates for long chain lengths. The decrease of the jam- 
ming density results from the formation of rigid semiloops 
in the granular chains which permit stable voids to exist 
in the packing, in excellent agreement with recent 3D ex- 
periments on granular polymers. For fixed chain length, 
the jamming density decreases when the chains are made 
stiffer since the rigid semiloops, and the voids stabilized 
by them, are larger. The jamming density can also be fur- 
ther decreased by increasing the fraction of granular loops 
present in the packing, which is also in agreement with 
experimental observations. The fact that our 2D simula- 
tions agree so well with the 3D experiments of Ref. 22 1 



indicates that the formation of semiloops in the chains 
is essentially a 2D phenomenon. In comparison to bidis- 
perse disk systems which show a nearly linear increase in 
the pressure as a function of density, characteristic of a 
critical phenomenon, in the granular polymer systems the 
pressure increases as a power law with exponent signifi- 
cantly larger than 1, suggesting that the jamming tran- 
sition in the granular chain system is different in nature 
from jamming in the bidisperse disks and may be related 
to the type of jamming that occurs for frictional grains. 
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